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ABSTRACT 
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We report evidence for cyclotron resonance features from the Soft Gamma Repeater SGR 1806 20 in 
outburst, detected with the Rossi X-ray Timing Explorer in the spectrum of a long, complex precursor 
that preceded a strong burst. The features consist of a narrow 5.0 keV absorption line with modulation 
near its second and third harmonics (at 11.2 keV and 17.5 keV respectively). The line features are 
transient and are detected in the harder part of the precursor. The 5.0 keV feature is strong, with an 
equivalent width of ~ 500 eV, and a narrow width of < 0.4 keV. Interpreting the features ^ electron 
cyclotron lines in the context of accretion models leads to a large mass-radius ratio (M/R > 0.3 M 0 /km) 
that is inconsistent with neutron stars, or requires a low (5-7) x 10 u G magnetic field that is unlikely for 
SGRs. The line widths are also narrow compared to those of electron cyclotron resonances observed so tar 
in X-ray pulsars. In the magnetar picture, the features are plausibly explained as ion cyclotron resonances 
in an ultra-strong magnetic field, which have recently been predicted from magnetar candidates. In this 
view, the 5.0 keV feature is consistent with a proton cyclotron fundamental whose energy and width are 
close to model predictions. The line energy would correspond to a surface magnetic field of 1.0 x 10 C 
for SGR 1806-20, in good agreement with that inferred from the spin-down measure in the source. 

Subject headings: Stars: Individual (SGR 1806-20) — Stars: Magnetic Fields - Stars: Neutron - 
X-Rays: Bursts 


1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) are a unique class of 
slowly rotating pulsars (P ~ 5 - 8 s) that glow quietly in 
X-rays (luminosity L ~ 10 35 — 10 36 ergs 1 ) for several 
years and suddenly become vigorously active for a few 
weeks to months, emitting hundreds of short (~ 0.1 s), 
bright {L ~ 10 39 - 10 42 ergs -1 ) bursts of soft 7-rays (see 
e.g. Hurley 2000 for reviews). Occasionally, SGRs also 
emit giant bursts that last for up to a few hundred sec- 
onds and exhibit remarkable pulsations that reveal their 
spin periods and confirm their nature as rotating neutron 
stars (Mazets et al. 1979; Hurley et al. 1999a; Ibrahim et 
al. 2001). 

The lack of evidence for a binary companion or a rem- 
nant accretion disk made the energy source of SGRs elu- 
sive. Two competing models offer contrasted views on 
SGRs as conventionally magnetized ( B « 10 12 G) neu- 
tron stars powered by fossil accretion disks (Marsden et 
al. 2001), or ultra-magnetized ( B ss 10 15 G) neutron stars 
powered by their own magnetic field (magnetar s; Duncan 
& Thompson 1992). Recently, strange/quark star models 
were also proposed for SGRs (Cheng & Dai 1998; Zhang, 
Xu & Qiao 2000). 

SGR 1806-20 is one of the four confirmed SGRs; three 
lie within the Galaxy (SGR 1900+14, SGR 1806-20, and 
SGR 1627-41) and one in the Large Magellanic Cloud 
(SGR 0526-66). The source was first detected by the Prog- 


noz 7, ISEE and KONUS instruments in 1979 (Laros et al. 
1986). A few years later, it underwent several periods of 
burst activity that were observed by different missions; 
this allowed a fairly accurate determination of the source 
position (Atteia et al. 1987; Laros et al. 1987). The lat- 
est episode of activity was closely monitored in 1996 with 
the Rossi X-ray Timing Explorer (RXTE). During these 
observations the source in quiescence was found to pul- 
sate with a 7.47 s spin period and to spin-down at a high 
rate (2.6 msyr -1 ). Interpreting the spin-down as due to 
magneto-rotational dipole losses leads to a magnetic field 
of ~ 8 x 10 14 G and a characteristic age of ~ 1500 yrs, typi- 
cal of a young pulsar (Kouveliotou et al. 1998) . SGR 1806- 
20 has been associated with the Galactic radio supernova 
remnant (SNR) G10.0-0.3 about 14.5 kpc away (Kulkarni 
et al. 1994; Corbel et al. 1997); however, this association, 
like those of other SGRs with SNRs (see Gaensler et al. 
2001), was later questioned on the basis of IPN and Chan- 
dra observations (Hurley et al. 1999b). Recently, a possible 
infrared counterpart was reported on the basis of Chandra 
observations (Eikenberry et al. 2001). 

To date, while a great deal of evidence has gathered in 
favor of the magnetar model over other scenarios, the de- 
termination of the magnetic field strength in SGRs (Kou- 
veliotou et al. 1998; Kouveliotou et al. 1999), although 
compelling, is still indirect (Marsden et al. 1999). Spectral 
signatures, on the other hand, promise decisive direct mea- 
surements of the field strength. In this Letter, we report 
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evidence for absorption features in the X-ray spectrum of 
SGR 1806 20. The nearly harmonic line spacing is sugges- 
tive of a cyclotron origin. We use the line energy to derive 
the star’s magnetic field and discuss the implications of 
this finding on current SGR models. 

2. OBSERVATION AND DATA ANALYSIS 

SGR 1806-20 entered an intense phase of bursting activ- 
ity in November 1996 during which it was extensively mon- 
itored with RXTE. Proportional Counter Array (PCA) 
data were extracted from the HEASARC archives and the 
high-resolution event mode data were used to construct 
the burst light curves. Among the bursts detected, we 
found one event that shows an unusual temporal profile 
with a long (« 0.5 s) multi-peak precursor followed by a 
bright burst (see Fig. la, b). The time history of this 
event shares a number of similarities with a remarkable 
burst from SGR 1900+14 where a 6.4 keV emission line 
was recently discovered in the burst precursor (Strohmayer 
& Ibrahim 2000; Ibrahim et al. 2001)= In fact: i) both 
bursts are preceded by a complex multi-peaked precursor 
and both precursors have similar temporal profiles that 
are different from typical SGR bursts; ii) both precursors 
are significantly longer than typical (~ 0.1 s) SGR bursts, 
iii) spectral evolution is detected in both precursors, while 
most bursts show a uniform spectrum. 

Motivated by these similarities we investigated in de- 
tail the spectral characteristics of the precursor. In order 
to subtract the background, we used a pre-burst persis- 
tent emission segment that is free from any bursting ac- 
tivity. We divided the precursor into four intervals and 
used the PCA data below 30 keV to fit the spectra of the 
different peaks separately. An absorbed power-law emis- 
sion model provided a qualitatively acceptable fit to the 
data and revealed a moderate spectral evolution as shown 
in Fig. lc. The spectrum is harder in the second inter- 
val with a larger y 2 than in the other three intervals. 
Absorbed bremsstrahlung models gave similar fits, with 
kT ~ 30 keV. The fit residuals (Fig. 2) show no signifi- 
cant structure in intervals 1, 3, and 4, but in interval 2, 
they suggest the presence of a narrow absorption feature 
near 5 keV. The fit was improved by adding a cyclotron 
absorption component, and the F-test gives a probability 
of 1.3 x 10 -3 for a chance reduction of the F-statistic by 
6.8 for a 5.0 keV absorption feature with a narrow width 
of 0.24 keV. The feature is significant at 3.2 a. 

Further inspection of the fit residuals suggested addi- 
tional modulation around the second and third harmonics 
of the 5.0 keV feature, and near 7 keV. The fit continued 
to improve with successive addition of three absorption 
lines. The energies derived from the fit are 7.5, 11.2, and 
17.5 keV respectively. With all four lines included, the F- 
test gives an improved random probability of 2.2 x 10 -4 for 
a chance reduction of the F-statistic by 6.2, corresponding 
to a 3.7 a significance. 

The spectrum and the predicted counts of the best-fit 
model are shown in Fig. 3 (right) along with the inci- 
dent photon spectrum that would be implied by the model 
(left). The fit parameters for the four intervals and the 
characteristics of the line fits are given in tables 1 and 2. 
The F-statistics given there are all in comparison to the 
absorbed power-law fit with no cyclotron resonance fea- 


tures. 

Absorption features did not improve the fits significantly 
for the other three intervals. We studied the implica- 
tions of instrumental effects in the PCA elsewhere and 
we showed that for count rates < 9 x 10 4 counts/s, pile- 
up and dead-time effects are not sufficient to modify the 
spectrum nor can they produce a spectral feature (Ibrahim 
et al. 2001; Strohmayer & Ibrahim 2000). The count rate 
during the main burst peaks exceeded this limit and we 
therefore excluded them despite some evidence for absorp- 
tion around 5 keV. The presence of the features in an in- 
terval with a moderate count rate and the fact that the 
features are present in part of the precursor and are not 
seen in intervals with comparable count rates and sim- 
ilar number of counts strongly argue against significant 
instrumental effects. Besides, when fitted with a gaussian, 
the 5.0 keV feature is strong, with an equivalent width of 
~ 500 eV, even larger than that of the 6.4 keV emission 
line from SGR 1900+14. This cannot be interpreted as 
due to imperfections in the PCA response matrix. Nu- 
merical simulations also showed that random fluctuations 
could produce power-law spectra with a y 2 as high as that 
of interval 2 in about 16% of the cases, but the fit residuals 
are randomly distributed and never showed a systematic 
structure similar to that of the real data. 

Taking into account that the features appear in one 
spectrum out of the four examined, the chance proba- 
bilities for the single feature and the set of features are 
5.2 x 10 -3 and 8.8 x 10 -4 , respectively. However, we have 
now looked at a large sample of bursts and found the 5 keV 
feature in the spectra of some other bursts as well. This 
will be reported in a follow-up paper (Ibrahim et al. 2002). 

Since an emission line has been detected in SGR 
1900+14 in a similar event, we checked this possibility also 
for present data. Visual inspection of the spectrum may be 
actually suggestive of a peak around 7 keV. However, the 
addition of a Gaussian emission line provided a marginal 
improvement in the fit with respect to an absorbed power- 
law, with a confidence level of only 74%. 

3. DISCUSSION 

The presence of electron cyclotron features in X-ray pul- 
sars was first predicted by Gnedin & Sunyaev (1974) and 
discovered a few years later by Triimper et al. (1978) and 
Wheaton et al. (1979). All observed cyclotron lines have 
been detected above 10 keV and are interpreted as electron 
features, with inferred magnetic fields B ~ (1-3) x 10 12 G 
(Cusumano et al. 1998; Heindl et al. 1999; Santangelo et al. 
1999). The presence of cyclotron features in SGRs would 
resolve a host of issues concerning the nature of these pe- 
culiar objects. SGRs have been considered as either ultra- 
magnetized neutron stars powered by the dissipation of 
their own B-field (Duncan & Thompson 1992) or conven- 
tional neutron stars powered by accretion from a very faint 
companion or a fossil disk (Marsden et al. 2001). Here we 
discuss the implications of our observation in light of these 
models. 

In a strong magnetic field and depending on the physi- 
cal conditions in the emitting/absorbing matter, cyclotron 
features from electrons, protons and other ions could in 
principle be observed. Possibly atomic absorption lines 
could also be present; for example, red-shifted iron at 



5 keV, but no model of generating such high flux and 
high temperature at a layer deeper than absorbing heavy 
atoms has been proposed. In a typical neutron star 
with B « 10 12 G, as invoked in accretion-powered mod- 
els, the fundamental electron cyclotron resonance lies at 
Ee = 11.6(1 + z) _1 (5/10 12 G) keV, where (1 + z) _1 = 
(1 - 2 GM/Rc*) 1/2 is the gravitational red-shift factor at 
the star surface. Proton and alpha particle resonances, on 
the other hand, would be undetectable in X-rays due to 
their very low energy ( E p<fi e ~ 2 - 5 eV). Interpreting the 
5.0 keV feature as an electron cyclotron resonance origi- 
nating close to the surface of a typical B ~ 10 12 - 10 13 G 
neutron star would require a very large gravitational red- 
shift (1 + z)~ l < 0.4 (or 2 > 1.5), with mass-radius ratio 
M/R > 0.3M©/km. Such values are inconsistent with 
neutron stars, for which M/R < 0.23M©/km as imposed 
by causality (Lattimer et al. 1990). Interestingly, these 
values are consistent with a more compact form of matter 
such as strange/quark stars; however, no predictions for 
spectral features from such objects have been made yet 
(Cheng k Dai 1998; Zhang, Xu k Qiao 2000). Fields be- 
low 4 x 10 11 G are ruled out since they imply 2 < 0. For 
the 5.0 keV feature to be an electron cyclotron resonance 
from a neutron star with acceptable mass and radius, the 
surface magnetic field could only be in the narrow range 
(5 _ 7) x 10 11 G, in apparent contrast with plausible val- 
ues for SGR 1806-20 (B ~ 10 13 G) within the accretion 
model (see Fig. 2 in Rothschild et al. 2000) and well below 
the average value in ordinary radio pulsars. The pulsar B- 
field distribution shows that ~ 20% of the population has 
B < 6 x 10 u G (e.g. Hartman et al. 1997; Regimbau k 
de Freitas Pacheco 2001). The rarity of SGRs as opposed 
to the larger population of pulsars with such a magnetic 
field also argues against this possibility. 

Moreover, electron cyclotron features in the spectra of 
accreting pulsars show large thermal broadening (A E « 
a few keV) due to the low electron mass in comparison 
with those of ions. Although the physical conditions in 
SGRs are not necessarily the same as in X-ray pulsars, 
broadening of the cyclotron line is to be expected in any 
case for typical electron energies ss 1 — 10 keV. For the 5.0 
keV fundamental the estimated width would be A E « 1.5 
keV. On the contrary, the observed feature is much nar- 
rower (AE « 0.24 keV). The 1 keV resolution of the PC A 
detector at 5 keV makes it difficult to precisely resolve fea- 
tures with widths < 1 keV. However, the 90% upper limit 
on the 5.0 keV feature width is 0.4 keV, still narrower than 
any known electron feature. All of the preceding consid- 
erations indicate that the interpretation of the observed 
5.0 keV feature in terms of an electron cyclotron line seems 
unlikely, although it can not be completely ruled out on 
the basis of present data. 

In the magnetar model ( B > 10 14 G), electron cy- 
clotron lines lie at very high energy E e > 1 MeV, out 
of range of the RXTE/PCA, whereas proton and alpha 
particle resonances come within reach, with fundamen- 
tals at E p = 6.3(1 + 2)- 1 (R/10 15 G) keV and E He = 
3.2(1 + 2)- 1 (B/10 15 G) keV. In this picture, the 5.0 keV 
feature is plausibly a proton cyclotron fundamental with 
an implied magnetic field B = 7.9 x 10 14 ( 1 + z) G. For a 
typical neutron star with M = 1.4 M© and R = 10 km, the 
surface field strength in SGR 1806-20 is B = 1.0 x 10 43 G, 
in a very good agreement with the value (B ~ 8 x 10 G) 


derived from the spin-down measure (kouveliotou et al. 
1998). The features at 11.2 keV and 17.5 keV are con- 
sistent with the second and third harmonics but with a 
slight deviation that could be due to emission from a re- 
gion with different magnetic field and/or red-shift. Elec- 
tron cyclotron features with slightly anharmonic spac- 
ing have been seen in transient X-ray sources (Heindl 
et al. 1999). Another possibility is that the 17.5 keV 
line is associated with the proton spin-flip transition at 
E$ = 17.6(1 + z) _4 (B/10 45 G) keV (Thompson 2000; 
Zane et al. 2001). However, this interpretation seems less 
plausible, since the rates for spin-flip transitions should 
contain an additional factor E p /m p c 2 1 with respect to 
transitions with conserved spin projection. 

In recent years the growing evidence for the existence of 
magnetars has stimulated several studies on the presence 
of ion cyclotron features in the X-ray spectra of SGRs 
(Bezchastnov et al. 1996; Thompson 2000; Zane et al. 
2001; Ho k Lai 2001). Our results on the energy and 
width of the fundamental feature are consistent with es- 
timates by Zane et al. (2001). However, the continuum 
underlying cyclotron features reported here is either non- 
thermal or bremsstrahlung for kT > 20 keV and probably 
originates in the star’s magnetosphere. In this respect, ion 
cyclotron lines may naturally form in the magnetar model 
of Thompson k Duncan (1995) as primary photons pro- 
duced in the pair fireball cross an optically thick baryon 
loaded sheath, located at < 27?. The presence of hard pho- 
tons with energy greater than the fundamental harmonic 
would also enhance the transitions between Landau lev- 
els and thus makes cyclotron features stronger (Gnedin k 
Sunyaev 1974). 

Similarly to proton cyclotron resonances, alpha parti- 
cle resonances could also be present. For the same value 
of B as implied by the observed proton fundamental, the 
alpha particle second harmonic would overlap the proton 
fundamental and the third harmonic would be consistent 
with the observed feature at 7.5 keV . The alpha particle 
fundamental would lie at 2.4 keV, too close to the lower 
end of the RXTE/PCA energy range to be separated from 
low-energy cut-off effects. Further observations with in- 
struments possessing higher spectral resolution and lower 
bandpass energy could resolve the alpha particle funda- 
mental and confirm or disprove this conjecture. 

Cyclotron features offer a new diagnostic tool for prob- 
ing SGRs. Further analysis and new observations promise 
to shed new insights and broaden our understanding of 
these enigmatic objects. This will potentially open new av- 
enues for testing quantum electrodynamics predictions and 
exploring new physical effects unique to ultra-strong mag- 
netic fields, which is one of the main reasons why SGRs 
are distinctively important. 
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Table 1 


Spectral Results For the Precursor 


Interval 

Photon Index 

N n (10 22 cm" 2 ) 

X 2 /do f 

F-statistic/Confidence level 

1 

1.8 ± 0.2 

21.4 ± 4.1 

35/32 


2 

1.4 ±0.1 

18.8 ± 3.7 

40/32 


3 

1.6 ±0.1 

18.2 ± 3.5 

12/32 


4 

1.9 ± 0.2 

26.7 ±7.5 

31/32 


2 (1 line) 

1.5 ± 0.2 

18.7 ± 3.4 

23.5/29 

6.8 / 99.87% 

2 (4 lines) (a) 

1.6 ±0.2 

21.1 ± 4.9 

13/24 

6.2 / 99.97% 


Note. ^ ^ : The line widths are held fixed at the best-fit values shown in table 2. 



Table 2 

Best-Fit Parameters for the Line Features 


Line Feature 

Energy (keV) 

Width (keV) 

Depth 

1 

5.0 ± 0.2 

0.24 ±0.1 

1.9 ±0.6 

2 

7.5 ±0.3 

0.45 ±0.2 

1.2 ±0.4 

3 

11.2 ± 0.4 

1.2 ±0.5 

0.9 ± 0.3 

4 

17.5 ± 0.5 

1.1 ±0.7 

1.0 ±0.4 


Note. — Cyclotron absorption model is described in Mihara, 
T. et al. 1990, Nature, 346, 251 
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FIG. 1.- Time history of the outburst from SGR 1806-20 as seen by RXTE/PCA (2-60 keyjThetoppanel (a)sho^ ; 
preceded by a long, complex precursor. Panel (b) shows the precursor intervals used m the spectral analys.s. Panel (c) sho 
of the power-law photon index for the precursor intervals. T 0 {UTC) - 9:13:43 on 1996 Nov. 18. 
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Fig. 2. The residuals of the four precursor intervals when fitted with an absorbed power-law model, 
interval 2 but are random in intervals 1, 3 and 4. 


Energy (keV) 

The residuals are significant in 




Energy (keV) 


Energy (keV) 


T* 16 spectrum and best-fit continuum model for the second precursor interval, with four absorption lines (RXTE/PCA 2-30 keV) 
(Right) The pulse height spectrum with the model predicted counts (histogram). (Left) The model (histogram) and the estimated photon 
spectrum tor the best ht model. 









